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X-RAY OBSERVATIONS OF PARSEC-SCALE TAILS BEHIND TWO MIDDLE-AGED PULSARS 
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ABSTRACT 

Chandra and XMM-Newton resolved extremely long tails behind two middle-aged pulsars, J1509-5850 and J1740+1000. The 
tail of PSR J1509-5850 is discernible up to 5.6' from the pulsar, which corresponds to the projected length l± = 6.5^4 pc, where 
d = Adi kpc is the distance to the pulsar. The observed tail flux is 2 x 10"^^ erg s"^ cm"^ in the 0.5-8 keV band. The tail 
spectrum fits an absorbed power-law (PL) model with the photon index F = 2.3 ± 0.2, corresponding to the 0.5-8 keV luminosity 
of 1 X IQ^^dl ergs s"', for = 2.1 x lO^^ cm'^. The tail of PSR J1740+1000 is firmly detected up to 5' ^ 2diA pc), with 
a flux of 6 X 10"''* ergs cm"^ s"' in the 0.4-10 keV band. The PL fit to the spectrum measured from a brighter, 3'-long, portion 
of the tail yields F = 1.4-1.5 and nn ~ 1 x 10^' cm~^; its 0.4-10 keV luminosity is ~ 2 x lO^^df^ ergs s~'. The luminosity of 
the entire tail is likely a factor of 3-4 higher. The large extent of the tails suggests that the bulk flow in the tails starts as mildly 
relativistic downstream of the termination shock, and then gradually decelerates. Within the observed extent of the J1509-5850 
tail, the average flow speed exceeds 5,000 km s"', and the equipartition magnetic field is a few x 10"^ G. For the J1740+1000 tail, 
the equipartition field is a factor of a few lower. The harder spectrum of the J1740+1000 tail impUes either less efficient cooling 
or a harder spectrum of injected electrons. For the high-latitude PSR J1740+1000, the orientation of the tail on the sky shows 
that the pulsar is moving toward the Galactic plane, which means that it was bom from a halo-star progenitor. The comparison 
between the J1509 and J1740 tails and the X-ray tails of other pulsars shows that the X-ray radiation efficiency correlates poorly 
with Ihe pulsar spin-down luminosity or age. The X-ray efficiencies of the ram-pressure confined pulsar wind nebulae (PWNe) 
are systematically higher than those of PWNe around slowly moving pulsars with similar spin-down parameters. 

Subject headings: pulsars: individual (PSR B 1509-5850, PSR J1740+1000) — stars: neutron — X-rays: stars 



L INTRODUCTION 

The multiwavelength synchrotron emission from ultrarel- 
ativistic particles produced in the pulsar magnetosphere and 
shocked in the ambient medium is observed as a pulsar-wind 
nebula (PWN; Rees & Gunn 1974; Kennel & Coroniti 1994; 
Arons 2007; Kirk & Lyubarsky 2007). Thanks to the high 
sensitivity and angular resolution of the Chandra and XMM- 
Newton observatories, about 50 PWNe have been detected 
in X-rays (Kaspi et al. 2006; Gaensler & Slane 2006; Kar- 
galtsev & Pavlov 2008, hereafter KP08). The X-ray obser- 
vations have shown that PWNe have complex morphologies, 
including toroidal structures around the pulsar, jets along the 
pulsar's spin axis, and cometary tails (see KP08 for Chan- 
dra images). In particular, if a pulsar moves with a super- 
sonic speed, V 3> Cj, the ram pressure, ;?ram = PambV^, exceeds 
the ambient medium pressure, Pamb^ resulting in a bowshock 
PWN with a tail behind the pulsar. The well-known exam- 
ple of such a bowshock-tail PWN is "the Mouse PWN" pro- 
duced by the young, energetic pulsar J1747-2958 (r = 26 kyr, 
£■ = 2.5 X 10"''' ergs s~'), which shows an X-ray tail of the pro- 
jected length Zj^ ^ 1 pc (Gaensler et al. 2004). In addition 
to the Mouse, about a dozen X-ray PWNe with bowshock- 
tail morphologies have been discovered recently (KP08). In- 
terestingly, one of the longest tails, l± ^ 1.5 pc, was found 
behind the 3 Myr old pulsar B 1929+ 10, with a relatively 
low spin-down power, £ = 3.8 x lO"*^ ergs s"' (Wang et al. 
1993; Becker et al. 2006; Misanovic et al. 2007). Although 
most of the known PWNe are powered by much younger and 
more powerful pulsars, it has been noticed (Kargaltsev et al. 
2007) that bowshock PWNe with extended tails are generally 
brighter than PWNe around slowly moving pulsars. 

The shape of the shock, the length of the tail, and the over- 
all appearance of the entire bowshock-tail PWN are expected 
to depend on the interplay between E, Pamb, Pamh, and v, as 
well as on the wind magnetization parameter a and the angle 



between the pulsar's spin axis and the direction of the pul- 
sar's motion. Analytical and numerical modeling of magne- 
tized outflows in bowshock-tail PWNe has been done by Ro- 
manova et al. (2005) and Bucciantini et al. (2005; B05 here- 
after), respectively. It was assumed in these works that the 
magnetic field in the tail is purely toroidal, and the pulsar's 
velocity is aligned with the rotation axis. In addition, the nu- 
merical calculations by B05 assumed a spherically symmet- 
ric relativistic pre-shock wind with an arbitrary magnetization 
parameter in the framework of ideal magnetohydrodynamics 
(MHD), while Romanova et al. (2005) postulated an equipar- 
tition between the magnetic and particle energy densities and 
included magnetic field reconnection into consideration. Be- 
cause of computational difficulties, the numerical calculations 
have been so far limited to short distances (a few shock stand- 
off radii) from the pulsar. To understand the large-scale prop- 
erties of extended tails and facilitate their modeling (including 
the intrinsic anisotropy of the wind and misalignment between 
the pulsar velocity and spin), observations of these objects are 
particularly important. 

Here we report on X-ray observations of very extended pul- 
sar tails behind two middle-aged pulsars, J1509-5850 and 
J1740+1000. PSR J1509-5850 (hereafter J1509) was discov- 
ered in the Parkes Multibeam Pulsar Survey' (Kramer et al. 
2003). The pulsar is located in the Galactic plane at the dis- 
persion measure (DM) distance of about 4 kpc (see Table 1 
for the pulsar properties). To study the X-ray properties of 
the pulsar and look for its PWN, we observed this field with 
Chandra. In addition to detecting the pulsar, the observation 
revealed a spectacular long tail, first reported by Kargaltsev et 
al. (2006). It was briefly described by Hui & Becker (2007; 
hereafter HB07), who also reported the discovery of diffuse 
radio emission, possibly associated with the tail. In this pa- 
per, we provide a detailed analysis of our Chandra observation 

' httpV/www.atnf.csiro.au/research/pulsar/pmsurv 
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TABLE 1 

Observed and derived pulsar parameters 



Parameter 


J1509-5850 


J 1740+ 1000 


R.A. (J2000) 


15''09'"27fl3(3) 


17''40"25?950(5) 


Decl. (J2000) 


-58°50'56f'l(5) 


+10°00'06^'3(2) 


Epoch of position (MJD) 


51,463 


51,662 


Galactic longitude 


319.97° 


34.01° 


Galactic latitude 


-0.62° 


20.27° 


Spin period, P (ms) 


88.9 


154.1 


Period derivative, P (10""*) 


0.92 


2.15 


Dispersion measure, DM (cm~' pc) 


137.7 


23.85 


Distance", d (kpc) 


3.8 


1.4 


Distance from the Galactic plane, z (kpc) 


0.04 


0.48 


Surface magnetic field, Bj (lO^^ G) 


0.91 


1.84 


Spin-down power, E (10^' erg s~') 


5.1 


2.3 


Spin-down age, t = P/(2P), (kyr) 


154 


114 



Note. — Based on the data from Kramer et al. (2003) and McLaughlin et al. (2002). Figures in parentheses represent la uncertainties in least-significant 
digits quoted. 



^ The distances are based on the dispersion measure and the Galactic electron density distribution model by Taylor & Cordes (1993). The Cordes & Lazio 
(2002) model gives d = 2.6 and 1.2 kpc for J1509 and J1740, respectively. 



of the tail and discuss implications of our results. 

Pulsar J1740-H1000 (hereafter J1740) has been discovered 
in an Arecibo survey (McLaughlin et al. 2000), and the 
follow-up radio observations were carried out by McLaugh- 
lin et al. (2002). The pulsar's properties are summarized in 
Table 1. At the DM distance of 1.4 kpc and Galactic lati- 
tude of 20.3°, J 1740 belongs to a very small group of pul- 
sars located at large distances from the Galactic plane, imply- 
ing either a progenitor from the halo population or a neutron 
star (NS) ejected from the Galactic disk with an exceptionally 
high speed, ~ 4000 km s~'. J 1740 was observed with Chan- 
dra in 2001 (PI: Z. Arzoumanian). In a preliminary analysis 
of these data we detected the pulsar, but the exposure was 
too short to study possible extended emission and the pulsar's 
spectrum. Therefore, we observed the field around J1740 in 
a much deeper imaging exposure with XMM-Newton. In this 
paper we report the discovery of the PWN associated with 
J1740 and focus on the analysis of the PWN properties, while 
the spectral and timing properties of the pulsar, including the 
detection of the thermal emission from the NS surface, will be 
presented in a separate publication (Misanovic et al., in prep.). 

Our X-ray observations of the J 1509 and J1740 tails and the 
data analysis are described in §2. We estimate the pulsar and 
flow velocities, the magnetic field strengths and energetics of 
the tails, compare these tails with the other pulsar tails cur- 
rently known, and discuss the implications of our findings in 
§3. Our main results are summarized in §4. 

2. OBSERVATIONS AND DATA ANALYSIS 

J 1509 was observed with the Advanced CCD Imaging 
Spectrometer (ACIS) on board Chandra on 2003 February 9 
(ObsID 3513). The useful scientific exposure time was 39.6 
ks. The observation was carried out in Very Faint mode, and 
the pulsar was imaged on the S3 chip, 35" from the optical 
axis (CHIPX=285, CHIPY=503). The other chips activated 
during this observation were SI, S2, S4, 12, and 13. The de- 
tector was operated in Full Frame mode, which provides time 
resolution of 3.24 s. For the analysis, we used the data re- 
processed on 2006 July 20 (ASCDS ver. 7.6.8, CALDB ver. 
3.2.4). 

J 1740 was observed with the Chandra ACIS on 2001 Au- 
gust 19 in Timed Exposure (imaging) mode (ObsID 1989; 



useful exposure time 5.1 ks) and Continuous Clocking mode 
(ObsID 2426; usefiil exposure time 16.5 ks). Only the former 
observation is useful for studying faint, extended features; it 
was carried out in Faint mode with the pulsar imaged on the 
S3 chip, 35" from the optical axis. The other activated chips 
were S2, 10, II, 12, and 13, and the time resolution was 3.24 
s. For the analysis, we used the data reprocessed on 2006 
December 21 (ASCDS ver. 7.6.9, CALDB ver. 3.2.4). 

J1740 was also observed with the European Photon Imag- 
ing Camera (EPIC) on board XMM-Newton on 2006 Septem- 
ber 28 and 30 (observations 0403570101 and 0403570201). 
The EPIC PN camera was operated in Small Window mode 
with the reduced field-of-view (FOV) of 4' x 4' and time- 
resolution of 6 ms, to analyze pulsations and perform phase- 
resolved spectroscopy of the pulsar. To obtain a large-scale 
image of the PWN, the EPIC MOS 1 and MOS2 cameras were 
operated in Full Window Mode, which provides a larger FOV 
at the expense of lower time resolution (2.6 s). The Thin filter 
was in front of the PN and MOS cameras during both observa- 
tions. The total exposures were 40.4 and 25.6 ks for the first 
and second observations, respectively. However, because of 
the reduced efficiency of the PN detector in the Small Window 
mode, the effective PN exposures were 28.6 ks and 17.8 ks in 
the first and second observations, respectively. We did not find 
significant flaring events in both observations (the background 
count rate varied by < 20%). Thus, the useful scientific expo- 
sures were 66 ks for MOS 1 and MOS2, and 46.4 ks for PN in 
the two observations combined. 

The Chandra data were reduced using the Chandra Inter- 
active Analysis of Observations (CIAO) software (ver. 3.4; 
CALDB ver. 3.4.0). The data obtained by XMM-Newton 
were reduced with the Scientific Analysis Software (SAS, ver. 
7.0). 

2.1. Images 

The X-ray images of the J 1509 and J 1740 fields reveal 
point sources at the radio pulsar positions. The best-fit cen- 
troid positions obtained from the Chandra ACIS images, are 
R.A. = 15''09'"27.163^ decl. = -58°50'56.12", and R.A. = 
17h4Qm25 939s^ decl. = +10°00'05.78", for J1509 and J1740, 
respectively. These positions differ from the corresponding 
radio positions by 0.3" and 0.6", respectively. The differences 



X-ray observations of parsec-scale pulsar tails 



3 





Fig. 1. — Large-scale ACIS view of the J1509's X-ray tail. The images show the S3 and S4 chips. The 1-7 keV energy band is chosen to maximize the S/N 
of the faint, extended tail. Both images are binned by a factor of 16 (pixel size is 7.8"). The image on the right is additionally smoothed with the 24" Gaussian 
kernel. "R" marks the counterpart of the pointlike source projected onto the putative radio tail (see §3.2 and Fig. 14). 



are comparable to the uncertainty of the absolute Chandra as- 
trometry^. 

The most striking features exposed in the binned Chandra 
and XMM-Newton images (see Figs. 1, 2, and 3) are long 
linear structures (we will call them "tails" hereafter) attached 
to J1509 and J1740. The tails have similar orientations on the 
sky, with the position angles of « 203° for J1509, and w 215° 
for J1740 (counted east of north). Assuming that the tails 
are caused by the pulsar motion, the corresponding position 
angles of the pulsars' proper motion can be estimated as w 23° 
and « 35°. 

2.1.1. J] 509 tail 

The ACIS image of J 1509 reveals the compact PWN (w 
10" X 30") in the vicinity of the pulsar (see Fig. 4). Because of 
its elongated, elliptical shape, we will call it the "head" here- 
after. Being the brightest close to the pulsar, the diffuse emis- 
sion gradually broadens and becomes fainter with increasing 
distance from the pulsar, forming an extended tail. The tail 
emission is firmly detected up to w 5.6' (or l± « 6.5ii4 pc, 
where d4 = (i/4kpc) (see Fig. 1) and possibly even farther, 
at a lower significance. This is demonstrated by the linear 
surface brightness profiles extracted along and across the tail 
(see Figs. 5 and 6). Although the collected number of counts 
from the tail is not very large (see Table 2), the low ACIS 
background still allows us to detect interesting morphological 
changes along the tail. For instance, beyond ^ 15" from the 

^ http://cxc. harvard. edu/caiydocs/cal_present_status.html#abs_spat_pos 



pulsar, the distribution of counts in the tail becomes asymmet- 
ric with respect to the inferred pulsar proper motion direction, 
shown by the dashed straight line in Figure 7 (e.g., there are 
more counts southeast of the line than northwest of the line). 
Being ^ 40" wide at the w 1.3' distance from the pulsar, the 
tail slightly brightens and then suddenly narrows by about a 
factor of two (as shown by the small arrows in the top panel of 
Fig. 7). Interestingly, very similar behavior is observed in the 
"Mushroom PWN" (KP08) associated with another middle- 
aged pulsar, B0355-I-50 (see Fig. 7, bottom, and the images in 
McGowan et al. 2006). Local tail brightenings are also seen 
at larger distances from the J1509 pulsar (Figs. 5 and 6); how- 
ever, the decreasing surface brightness of the tail and reduced 
off-axis resolution do not allow one to check whether or not 
these brightenings are accompanied by significant morpho- 
logical changes. 

2.1.2. J] 740 tail 

Even in the short, 5 ks, ACIS imaging exposure of J 1740 
we can see traces of faint extended emission southwest of the 
pulsar in the heavily binned ACIS-S image (Fig. 8), and the 
unbinned image of the immediate vicinity of the pulsar shows 
an excess of counts in the same direction and perhaps some 
emission in front of the pulsar (see the inset in Fig. 8). The 
long (~ 5.5', or l± ^ 2c/i 4pc, where d-14 = d/lAkpc) tail is 
much better seen in the EPIC images (Figs. 2 and 3). Similar 
to J1509, the Hnear profile of the J1740 tail (see Fig. 9) hints 
at surface brightness fluctuations along the tail; however, de- 
tailed investigation is not possible because of the high EPIC 
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TABLE 2 

Count statistics for PWN regions 



Region 


A 


Mot 






S/N 


S 


J 1509 


Head 


0.073 


118±11 


18±1 


100± 11 


9.1 


34.6±3.8 


Rl 


1.16 


542 ±23 


249 ± 11 


293 ±26 


11.3 


6.38 ±0.57 


R2 


1.26 


453 ±21 


280± 12 


173 ±24 


7.2 


3.45 ±0.48 


R3 


2.24 


691 ±26 


530± 16 


161±31 


5.2 


1. 80 ±0.35 


Entire tail 


4.68 


1686±41 


1059 ±23 


627 ±47 


13.3 


3.38 ±0.25 


J1740 


Bright portion (MOSl+2) 


2.10 


1658±41 


1026 ±32 


632 ±56 


11.4 


6.48 ±0.57 


Bright portion (PN) 


2.10 


2445 ±49 


1791 ±42 


654 ±70 


9.2 


6.70 ±0.72 



Note. — Total (Mot), scaled background (Mg), and source (Mrc) counts are in the 0.5-8 keV band for the Chandra data, and 0.4-10 keV band for XMM- 
Newton combined data, extracted from the regions of area A (in arcmin") shown in Figs. 10 and 3 for J1509 and J1740, respectively. The mean surface brightness, 



cS, is in units of counts ks 



-- in the 0.5-8 keV and 0.4-10 keV bands for J1509 and J1740, respectively. 




Fig. 2.— Combined MOS1+MOS2 images (18' X 13'; 0.3-lOkeV) of the 
field around J1740 from two XMM-Newton observations (see text). The top 
panel shows a binned image (pixel size 8") while the bottom image (pixel 
size 4") is adaptively smoothed with a Gaussian kernel. The top panel shows 
the regions from which the linear brightness profiles (shown in Fig. 9) were 
extracted. The overlaid contours in the bottom image show optical and ra- 
dio sources from the DSS2-blue (shown in cyan), DSS2-red (magenta), and 
NVSS 1.4 GHz (green) images. 

background and poor angular resolution of XMM-Newton. In 
addition to the long tail, the smoothed MOS1+MOS2 image 
of the J 1740 field (Fig. 2) and the linear profile in the east- 
west direction (Fig. 9) show some diffuse emission w 2.5' east 
of the pulsar, whose connection to the pulsar remains unclear. 
Since the XMM-Newton resolution is much more coarse than 




Fig. 3. — 9.4' X 6.2' binned image (pixel size 8") of the field around J1740 
combining MOS1+MOS2 (Full Frame mode) and PN (Small Window mode) 
data of two XMM-Newton observations, in the 0.3— lOkeV band (see text). 
The regions used for spectral extraction are marked by an ellipse (the source 
region) and a circle (the background). Galactic coordinate grid is shown. 



that of Chandra, emission from a number of clustered point 
sources may sometimes mimic the diffuse emission in the 
EPIC images. Therefore, we searched for point sources at 
radio and optical wavelengths. The bottom panel of Figure 
2 shows the contours of optical and radio sources detected in 
the DSS2-blue, DSS2-red and NVSS (1.4 GHz) images of the 
region around J 1740. There are a few optical sources that co- 
incide with the X-ray pointlike sources in the EPIC images. 
However, the majority of optical sources in the field are faint 
and have no X-ray counterparts. Furthermore, the distribution 
of these sources shows that they cannot be responsible for the 
emission detected in the tail region, nor in the region « 2.5' 
east of the pulsar. Yet, until a high-resolution, deep ACIS im- 
age is obtained, we cannot exclude the possibility that some of 
the J1740 tail emission comes from faint point sources lacking 
obvious optical counterparts. 

As the orientation of the J1740 tail on the sky indicates 
the proper motion direction, it provides useful information 
about the pulsar's birthplace. Being located high above the 
Galactic plane, at z = 0.48iii.4 kpc, the pulsar might be a 
very high-speed NS if it were born in the plane (e.g., v > 
4200c/i 4(114kyr/f) km s"', where t is the pulsar age, if it 
were born in the midplane of the Galaxy). However, the X- 
ray image of the tail, shown together with the Galactic coor- 
dinate grid in Figure 3, suggests that the pulsar is moving at a 
small angle, ~ 7°, toward the plane, which means that it was 
born well out of the Galactic plane, likely from a halo-star 
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Fig. 4.— Top: ACIS-S3 image of J1509 and its vicinity (1-7 keV; pixel 
size 0"49). Also shown are the extraction regions used to measure the pulsar 
spectrum (r = 0.9" circle) and the spectrum of the PWN head (14" X 24" 
ellipse). The radio pulsar position is shown by the small black ellipse, whose 
size shows the position uncertainty. Bottom: Adaptively smoothed sub-pixel 
resolution image (pixel size 0.25") of the same region obtained by remov- 
ing the pipeline pixel randomization and subsequently applying the sub-pixel 
resolution tool (based on analyzing the charge distribution produced by an 
X-ray event; Tsunemi et al. 2001; Mori et al. 2001). The line shows the ap- 
proximate symmetry axis of the PWN, which can be crudely identified with 
the direction of pulsar's proper motion. 

progenitor. 

2.2. PWN spectra 
2.2.1. J 1509 

We extracted the spectrum of the J 1509 PWN from the three 
numbered rectangular regions shown in Figure 10. Region 
1 (Rl hereafter) includes the PWN head, but not the pulsar 
(we excluded counts from the 0.9" radius circle shown in Fig. 
4). We also extracted the spectrum from the combined region 
1+2+3 (although it includes the PWN head, the head contribu- 
tion is small compared to that of the tail, so that we will refer 
to it as the "entire tail" region). For each of the regions, the 
background was extracted from the two corresponding rectan- 
gular regions located on both sides of the source region (see 
Fig. 10). The measured source and background count num- 
bers and areas are given in Table 2. The observed flux in the 
entire tail is Ftaii = (2.0 ±0.2) x 10"'^ ergs s"' cm"^ in the 0.5- 
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TABLE 3 

PL FITS TO THE J 1509 AND J 1740 PWN SPECTRA 



Model 


"H,22'* 




r 


X^/dof 


^-X.32'^ 




6-8= 


J 1509 


Head 


[1.3] 


0.7 ±0.2 


1.27±0.25 7.56/9 


1.15±0.32 


23 


2.7 


Head 


[2.1] 


1 4-0.17 


1.83±0.30 


8.37/9 


1.33±0.41 


27 


5.3 


Rl 


[2.1] 


5 7+' 2 
-'•'-1.0 




3.60/8 


4.4 ±0.4 


5.5 


1.4 


R2 


[2.1] 


4 3+1.7 

^•-'-1.2 


2.36 ±0.40 


1.68/4 


2.9 ±0.6 


3.7 


0.95 


R3 


[2.1] 


4.7±1.6 


T ^^-0.43 
^■-"+0.35 


2.16/4 


3.2±0.7 


2.0 


0.58 


Entire tail 


[2.1] 


16.5 ±2.4 


2 36+'' l5 
^■^°-0.17 


9.1/14 


11.2± 1.0 


3.5 


1.0 


J 1740 


Bright tail 
Bright tail 


05+0 " 
"■"-'-0.05 

[0.1] 


n 73+0.25 
"■'^-0.16 
0.80±0.11 


1 ^^+o.20 

-0.25 

1.48 ±0.19 


50.9/38 
51.5/39 


0.16±0.04 
0.15±0.04 


0.91 
0.88 


0.1 
0.1 



Note. — The uncertainties are given at the 68% confidence level for one 
interesting parameter 

"The fits for J1509 are for fixed «h.22 = "h/IO^^ cm"^. For J1740, the fits 
are for both free and fixed nH.22 (fixed values are in brackets). 

''Spectral flux in units of 10"^ photons cm"- s"' keV"' at 1 keV. 
Isotropic luminosity in units of 10^- ergs s"' in the 0.5-8 keV band for 
J1509 and 0.4^10 keV band for J1740. 

''Mean unabsorbed intensity in units of Ur" ergs s"' cm"^ arcsec"^ in the 
0.5-8 keV band for the J 1509 tail, and in the 0.4-10 keV band for the J 1740 
tail. 

B-8 = B/( 10"* photons s"' cm"^ arcsec"^); B = M/A is the average spec- 
tral surface brightness at i? = 1 keV. 



8 keV band. In addition, we extracted the head spectrum from 
the elliptical region shown in Figure 4 (middle panel), with the 
absorbed flux Fhead = (4.9 ±0.5) x lO"'"* ergs s"' cm"^ in the 
0.5-8 keV band. 

First, we fit the Rl spectrum with the absorbed power-law 
(PL) model, allowing the hydrogen column density, nu, to 
vary. The fit gives 11^22 = ne/lO^^cm"^ = 2.1^q5 (the un- 
certainties here and below correspond to the 68% confidence 
level for a single interesting parameter, i.e., x^~Xmin - 
Fitting the same model to the R2 and R3 spectra results in 
larger uncertainties, while the best-fit «h values are similar to 
that obtained for Rl. The PL fit to the spectrum of the entire 
tail gives a slightly larger column density, nH.22 = 2.5 ± 0.5, 
while the fit to the head spectrum gives a somewhat lower 
«H,22 = 1.3 ±0.5, which is still consistent with the above val- 
ues (see Fig. 11)^. The total Galactic HI column in that direc- 
tion is nHi.22 ~ 1.7 (Dickey & Lockman 1990); however, this 
does not contradict the larger ne values for the J 1509 PWN 
because the «h deduced from an X-ray spectrum under the as- 
sumption of standard element abundances generally exceeds 
the hhi measured from 21 cm observations by a factor of 1.5- 
3 (e.g., Baumgartner & Mushotzky 2006). Finally, the value 
nH,22 = 0.42, estimated from the pulsar's dispersion measure 
under the standard assumption of 10% ionization of the in- 
terstellar medium (ISM), is noticeably lower than the above 
values, a which possibly means that average ISM ionization 
is lower than 10% in the direction to J 1509. Thus, given dif- 
ferent kinds of uncertainties, we can only conclude that the 
true «H,22 is somewhere between 0.4 and 3. Throughout most 
of the analysis below we choose to fix «h,22 at the value ob- 
tained from the fit to the Rl spectrum since it has the highest 

' These values are larger than nH.22 = 0.82+jj reported by HB07 from the 
same ACIS data. From the description in HB07, it appears that the lower hh 
was a result of overcon'ecting for the ACIS filter contamination. To correct 
for the contamination, HB07 used the XSPEC model ACISABS despite the 
fact that this correction is already included in the ARF files generated by 
CIAO (since the release of CALDB 2.26 on 2004 February 2). 
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Fig. 5. — ACIS images (top) of the J1509 tail (binned to the pixel size of 4" and smoothed with a 20" Gaussian kernel) in the total (1-7 keV; lefi), soft (1-2.5 
keV; middle), and hard (2.5-7 keV; right) bands. The total band image is color coded, with the soft emission shown in red and hard emission shown in green. 
The bottom panels show the tail profiles in each energy band. The linear brightness profiles are obtained by extracting the source and background counts from 
the 1.6' X 6' rectangular regions shown in the top left panel. The original unsmoothed images were used for the extraction. The profile extracted from the source 
region is plotted by the solid line while the averaged background profile extracted from the adjacent background regions is shown by the dashed line. The bin 
size is 15" in the linear profile plots. 



S/N among the other PWN elements (see Table 2)^. 

With the fixed n^^i = 21, the absorbed PL fits to the Rl 
spectrum and the entire tail spectrum result in the photon in- 
dices Tri = 2.1 ±0.2 and Ttau = 2.4 ±0,15, with the unab- 
sorbed 0.5-8 keV fluxes F™"''^ = 2.33!j;:g x 10"'^ ergs cm'^ 
s-i and F™*' = (5.9±0.5) x 10"'^ ergs cm'^ s"' (see also Fig. 
12 and Table 3). The unabsorbed 0.5-8 keV luminosity of the 
entire tail is Ltaii ~ 1.1 x \(fi^dl ergs s"', while the luminos- 
ity of the head is a factor of 10 lower (see Table 3 for other 
details). 

The pulsar's spectrum will be analyzed in detail elsewhere. 
For the purposes of this paper, it suffices to say that the spec- 
trum fits the absorbed PL model with Fpsr — 2.2 and an un- 



absorbed 0.5-8 keV flux F^f'' 



(for «H,22 



7.8 X 10 ''^ ergs s ' cm ^ 
2.1), corresponding to the isotropic luminosity 



Lpsr=l-5x 10^2^1 ergs s"'. 

2.2.2. 77740 

Since the ACIS imaging exposure of J 1740 was very short, 
we have to rely upon the XMM-Newton data to infer the spec- 
tral properties of the extended tail. We excluded events with 
energies below 0.4 keV and above 10 keV from the analy- 
sis because of the large particle background contribution at 
these energies. We then extracted the spectrum of the bright 
part of the tail (within the PN FOV) from the elliptical re- 
gion shown in Figure 3. The choice of the background region 
was limited by the size of the PN FOV in the Small Window 

* Although the spectrum extracted from the entire tail region has even 
higher S/N and formally fits with a PL model, the actual spectrum of the entire 
tail can deviate from a simple PL because of synchrotron cooling. (Evidence 
of synchrotron cooling is seen in the energy-resolved images shown in Fig. 
5.) 



mode, so we extracted it from the largest possible circular re- 
gion, also shown in Figure 3, In the first observation we col- 
lected 1459, 471, and 492 counts (of which 25,6%, 38,0%, 
and 27,9% counts come from the source) in PN, MOSl, and 
MOS2, respectively, A slightly smaller numbers of counts 
were collected from the second observation: 986 (28,4% from 
the source) in PN, 342 (46,1%) in MOSl, and 353 (44,9%) in 
MOS2, 

We fit the PN and MOS spectra of both observations si- 
multaneously in the 0,4-10 keV band using an absorbed PL 
model. First, we fit the data with all parameters allowed to 
vary. The hydrogen column density obtained from this fit is 
poorly constrained, «h.22 = 0-0,16 (see Fig, 13, top, and Ta- 
ble 3), Then, we fixed the absorption column to the value 
nH,22 = 0,1, consistent with the pulsar's dispersion measure 
assuming the 10% ISM ionization and the fits of the pul- 
sar spectrum (Misanovic et al,, in prep,). The results of the 
spectral fitting are shown in Table 3 and Figure 13 {bottom). 
The energy flux from the tail is Ftaii = (6,0 ±0.6) x lO"''* ergs 
cm"^ s"' in the 0.4-10 keV band. The tail spectrum is rel- 
atively hard, with the photon index V^mX = 1-1-1.7; its unab- 
sorbed flux, /^ta™*"* ~ 6.6 X 10"'"* ergs s"' cm"^, corresponds 
to the 0.4-10 keV luminosity Ltaii w 1.5 x lO^^dj^ ergs s"' 
(Ltaii = 1.3 X IQ^^dl^ergs s"' in the 0.5-8 keV band). 

The MOS and ACIS images (Figs. 2 and 8) suggest that 
the tail extends beyond the elliptical region used in the above 
spectral analysis, up to ^ 7' from the pulsar. However, poor 
statistics precludes meaningful analysis of the spectra ex- 
tracted outside the PN FOV A crude estimate of the remaining 
tail flux can be obtained from net tail counts measured in the 
MOS and ACIS images outside the elliptical region. Using 
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Fig. 6. — Top: ACIS-S3 image (1-7 keV) showing tlie briglitest part of the 
J1509 tail. Bottom: Linear profiles across the tail obtained by extracting the 
source counts from the two rectangular (132" X 15") regions shown in the 
top image. The original unsmoothed 1-7 keV image was used for the extrac- 
tion. The profiles extracted from the regions closer to and farther from the 
pulsar are plotted by the black and red histograms, respectively. The dashed 
black histogram shows the averaged background profile extracted from the 
two background regions on both sides of the tail. 

PIMMS^ and assuming «h.22 = 0.1 and the PL slopes inferred 
above, we estimated the luminosity of the entire visible tail to 
be Ltaii ~ (2-5) X IQr'^dl^ ergs s"'. 

The detailed phase-resolved analysis of the pulsar spectrum 
shows that it includes both thermal and non-thermal emission 
components (the results will be reported in a subsequent pa- 
per; Misanovic et al., in prep.). Here we just note that the non- 
thermal part of the spectrum fits a PL model with Fpsr ~ 1 .2- 
1.4 and luminosity Lpsi- ~ 3 x lQ^^d\^ ergs s"' in the 0.5-8 
keV band. 

3. DISCUSSION 

The very long tails of the J 1509 and J1740 pulsars are ex- 
treme examples of ram-pressure confined outflows from su- 
personically moving pulsars. To date, about 13 extended X- 
ray sources likely associated with fast-moving pulsars have 
been found, although in some cases the interpretation of the 
observed elongated structures as tails is tentative (see KP08 
for a review). The sample of well-established X-ray PWNe 
with long ram-pressure confined tails is still very limited (see 
§ 3.3 and Table 5). We expanded the sample by adding the 
two pulsars, J1740 and J1509, with remarkably long, promi- 
nent X-ray tails. 

Despite the similarities in their appearance, pulsar tails 
show significant differences in their X-ray spectra, luminosi- 
ties, apparent lengths, surface brightness distributions, and 

^ http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 
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Fig. 7. — Top: Images of the J1509 tail (in the 1-7 keV band, binned to 
the pixel size of 2"). The parabola, x = 2.45"[(z + 0.5")/l"]'/-, where z is 
the distance from the pulsar along the tail, shows an approximate boundary 
of the PWN head within ~ 15" from the pulsar, but it does not fit the tail's 
shape at larger distances. The part of the other parabola, .r = 0.44" (z/ 1 " ) ' 
(shown by the dashed line), corresponds to eqs. (32-34) in Romanova et al. 
(2005), for v„,, = 10,000 km s"', Vpsr = 300 km s"' and fo = 0.5" (see also 
§3.2). The white arrows show an approximate observed width of the tail (see 
text). Bottom: The Mushroom PWN powered by PSR B0355+54. Notice the 
sudden narrowing (transition from the "cap" to the "stem"), similar to that 
seen in the J 1509 tail. 

multiwavelength properties. Below we derive physical prop- 
erties of the J 1509 and J1740 tails, compare them with model 
predictions, and discuss the physical processes that determine 
the tails' appearance and spectra. We also provide a com- 
parison with other pulsar tails, stressing the differences and 
discussing their possible origin. 

3.1. Morphology of the tails and pulsar speeds. 

The appearance of the X-ray PWNe around J1509 and 
J 1740 clearly shows that the pulsars move with supersonic 
speeds. If the pulsar speed is supersonic, the shocked pul- 
sar wind (PW) ahead of the pulsar is confined between the 




Fig. 8. — ACIS images of the J1740 field. Tlie main panel shows the 
binned (pixel size 7.9") ACIS image in the 0.5-6 keV band smoothed with 
a 24" Gaussian. Faint emission from the extended tail is visible within the 
region between the dashed lines. The inset shows the immediate vicinity 
of J 1740 in the same energy band. The image resolution is maximized by 
applying a sub-pixel resolution tool (Tsunemi et al. 2001; Mori et al. 2001). 
The pixel size in the inset image is 0.25". 
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Fig. 9. — Top: Linear profiles of the J1740 tail and the background obtained 
by extracting counts from the rectangular regions along the tail and a parallel 
source-free region located north-east from the tail (the two 1 ' X 11 .5' rectan- 
gular regions are shown in Fig. 2). Bottom: Linear profiles extracted along 
the east-west direction from the two 1.3' X 10.3' rectangular regions shown 
in Fig. 2. The unsmoothed MOSl-^MOS2 image in the 0.4-10 keV band, 
combining the two XMM-Newton observations, was used for the extraction. 



J 1509 




Fig. 1 1 . — Confidence contours (68% and 90%) in the nn-F plane for the 
PL fit to the J1509 head (blue), region Rl (green), and the entire tail (black) 
spectra. The contours are obtained with the PL normalization fitted at each 
point of the grid. The T values obtained from the fits with fixed «h.22 = 0-42 
(inferred from the DM value assuming 10% ISM ionization) are shown by 
the error bars. 
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Fig. 12.— Confidence contours (68% and 90%) for the PL fits to the 
J1509 PWN regions (labeled in the plot) for the fixed nH.22 = 2.1. The PL 
normalization is in units of 10"^ photons cm"- s"' keV"' at 1 keV. The dashed 
curves are the loci of constant unabsorbed flux in the 0.5-8 keV band; the flux 
values near the curves are in units of 10"'^ ergs cm"^ s"' . 
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Fig. 13.— Top: Confidence contours (68%, 90%, and 99%) in the hh- 
r plane for the PL fit to the J 1740 tail spectrum. The contours are obtained 
with the PL normalization fitted at each point of the grid. Bottom: 68%, 90%, 
and 99% confidence contours for the fixed nn = Ix 10^' cm"^. The dashed 
lines are the loci of constant unabsorbed flux in the 0.4-10 l£eV band; the flux 
values are in upUs of 10"'5 ergs cm"^,s^nr.o\ j »i, » » j- 

Dow-snapea terrmnatioii sliock (TS) and the contact disconti- 
nuity (CD) surface, while the shocked ambient matter is con- 
fined between the CD surface and the forward shock (FS). The 
PWN behind the pulsar acquires an elongated shape stretched 
in the direction opposite to the direction of the pulsar's mo- 
tion. Recent numerical simulations by BOS show that for the 
idealized case of an isotropic PW^ the TS acquires a buUet- 
like shape, with the TS apex (bullet head) at the distance 



[47rc(/?an,b + /5ram)J 

from the pulsar, where the PW pressure, p„ « £(47rcr^)"', 
is balanced by the sum of the ambient pressure, pamb, and the 
ram pressure, p^^m = pv^ = 1 .67 x lO'^'^nvj dyn cm"^ (here vv = 
v/lO^cms"' is the pulsar speed, and n = p/niu is in units of 
cm~^). Assuming pram ^ Pamb [equivalent to 1, where 

M. = v/Cj = (3pram/5pamb)''^^ is the Mach number] we obtain 

Rh « 4.0 X 10^^£35^«~'''^^7' cm, (2) 
or, equivalently, /^ram ~ 2.65 x 10~'^ EisR^^^f^ dyn cm~-^, where 
£35 = £/(1035 grgg g-i^^ ^ R^/Wcm. 

^ If the PW is mostly confined in the equatorial plane (i.e., in the plane 
perpendicular to the pulsar's spin axis), as we see in young pulsars, then the 
TS shape may strongly depend on the angle between the pulsar's velocity and 
the spin axis. 



For large Mach numbers and small values of the magnetiza- 
tion parameter a of the pre-shock pulsar wind, the TS bullet's 
cylindrical radius is rxs ~ Rh , and the distance of its back sur- 
face from the pulsar is Rb ^ 6Rh (BOS). The numerical sim- 
ulations (BOS) and analytical models (Romanova et al. 200S) 
also show that an extended tail forms behind the back surface 
of the bullet^. According to BOS, the flow speed within this 
tail ranges from (0.8-0.9)c in its sheath to (0.1-0.3)c in the 
narrow inner channel of the tail. 

For J1S09, the 3"-long region of enhanced surface bright- 
ness immediately behind the pulsar (see Fig. 4, top) can be 
interpreted as emission from the TS bullet, implying the dis- 
tance R/j ^ Rh/(3 ^ 3 X 10'*'((i4/ sini) cm to the (unresolved) 
TS apex (; is the angle between the pulsar velocity and the 
Une of sight). This assumption and equation (2) give ~ 
1 .S X 10"^(c/4 / sin dyn cm"^ and an estimate for the pulsar 
speed as a function of ambient number density, 

V w 400eI^\-^/^RI\^ km s"^ ~ 300n-^/^(J4/ sin «)"' km s"^ . 

(3) 

Alternatively, we can estimate the Mach number as a function 
of ambient pressure, 

M ^ (3£/207r/?^cpamb)'/' - mdAl sin/)-' p:^^^_^^, (4) 

where Pamb,-i2 = Pwi,/ 10"'^ dyncm"^ (Pamb,-i2 ~ 1 is a typical 
ISM pressure in the Galactic plane; e.g., Ferriere 2001). 

Unfortunately, the poor angular resolution of XMM- 
Newton and the shortness of the Chandra exposure do not al- 
low us to constrain the size of the TS bullet around J 1740. 
However, the Chandra image of the vicinity of J 1740 (see the 
inset in Fig. 8) shows a brightening ~ 0.7" ahead of the mov- 
ing pulsar One can speculate that the brightening corresponds 
to the TS apex, which gives Rh^ l.S x 10'^(c/i.4/ sin;) cm. 
This corresponds to the pulsar speed 

v~400n"^/2(i/i.4/sinO"^kms"^ (S) 
and the Mach number 

130(^1.4/ sin/)- V:iK-,3 (6) 
(notice that we changed the scaling here to account for the 
lower pressure at large distances from the Galactic plane). 
The estimated pulsar speed is consistent with v±_ ~ 200 km 
s"\ inferred from the interstellar scintillation measurements 
by McLaughlin et al. (2002), if « ^ Ad'^\(?,mi)'^ cm"^. It can 
be reconciled with the space-averaged density n ~ 0.04-O.OS 
cm-^ expected at z ~ 400-SOO pc (Fig. 2 in Ferriere 2001), if 

1/2 —1/2 

sin / ^ 0.3t/[ 4 , which would correspond to v ^ 600c/[ 4 km 

s-i, M ~ "^dp'^l-nd-'iT ' the tail's length / > 6d\']; pc. 

As the ISM pressure is much more uniform than the den- 
sity, the Mach number values (eqs. [4] and [6]) are more cer- 
tain than the speed values in equations (3) and (5). On the 
other hand, the ambient pressure in the pulsar vicinity may 
be higher than the typical ISM pressure because the pulsar's 
UV emission heats and ionizes the surrounding medium (e.g., 
Bucciantini & Bandiera 2001) . Therefore, the Mach numbers 
can be somewhat lower than the numerical values in equations 
(4) and (6). Also, we should not forget that our estimates are 
based on the assumption of isotropic distribution of the un- 
shocked pulsar wind. The estimates may change considerably 
if the wind is essentially anisotropic, but no models have been 
published for this case so far 

' Note that the numerical solutions by BOS do not extend beyond ~ lOi?/, 
from the pulsar. 
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Note. — Parameters used to estimate the equipartition magnetic field, 
Beq = Beq.-e from the measured surface brightness, = 0.1 

keV and Eu = 10 keV (see Table 3 and §3.2). The numerical coefficients 
ymp, yup > and Up are estimated using Table II in Ginzburg & Syrovatskii 
(1965). ' 

3.2. Magnetic field and fiow speed 
Assuming some value for the ratio km = Wmag/wrei of the 
magnetic energy density, Wmag = B^/(87r), to the energy den- 
sity of relativistic particles, Wrei, the magnetic field in the 
tails can be estimated from the measured synchrotron surface 
brightness (Pavlov et al. 2003): 

(7) 

Here B = Af/A =10 ^B-g photons (s cm^ keV arcsec^) ^ is the 
average spectral surface brightness at £ = 1 keV (see Table 
3), J\f is the normalization of photon spectral flux measured in 
area A, s= 10'^ sn cm is the average length of the radiating re- 
gion along the fine of sight, £„,,/> = E,n/ym,p, Em,p = -Eivi/yM.p, 
Ejji and Em are the lower and upper energies of the pho- 
ton power-law spectrum (in keV), and p, yM,p and Up are 
the numerical coefficients whose values depend on the slope 
p = 2T-l of the electron power-law spectral energy distribu- 
tion (see Table 4 and Ginzburg & Syrovatskii 1965). 

The quantities B and F have been directiy measured (Ta- 
ble 3), and s can be estimated assuming that it is close to 
the observed tail's width. In addition to these parameters, 
the strength of the magnetic field depends on the boundary 
energies of the photon power-law spectrum, which are rather 
uncertain because of the lack of sensitive high-resolution ob- 
servations outside the X-ray range. However, the observed 
softening of the J 1509 tail spectrum with increasing distance 
from the pulsar (Table 3), as well as the decrease of the tail's 
length at higher photon energies (Fig. 5), suggest that, at least 
for the J1509 tail, Em is close to the upper energy of the ACIS 
band, ~ 10 keV. The lower energy of the synchrotron spec- 
trum is not constrained; we can only be sure that is smaller 
than the lower energy of the ACIS band, ~ 0.3 keV. In the nu- 
merical estimates (Table 4), we will assume = 0.1 keV and 
Em = 10 keV. 

Under these assumptions, we obtain the equipartition mag- 

—2/7 

netic fields Beq = Bk,„ ^ 10-20 jiG in several regions of the 
J 1509 tail, with a hint of a decrease of fieq with increasing 
distance from the pulsar (see Table 4). For the bright portion 
of J 1740 tail, where the spectrum was measured, we found 

Beq ~ 5 fiG. 

We stress that the magnetic fields in Table 4 may differ by a 
factor of a few from the actual magnetic fields. First of all, if 
the energy range of the power-law spectrum is much broader 
than the 0. 1-10 keV band, the magnetic field is larger than the 
estimates in Table 4. For instance, as F > 1 .5 (p > 2) for the 
spectra of the J1509 tail, the derived magnetic field is sensitive 



to the unknown E^ value: B cx eJ ' . If the radio emis- 
sion in the SUMSS** image reported by HB07 (see also Fig. 
14) indeed belongs to the J 1509 tail, then the X-ray power- 
law spectrum might be extrapolated down to 10-4-10-2 
keV (depending on the F value), which would correspond to a 
factor of 1.5-6 higher magnetic fields. On the other hand, the 
magnetic field strength may be lower than Beq if km < 1 (e.g., 
by a factor of 1.9 if k„ = 0.1.) Overall, it seems reasonable to 
conclude that the magnetic fields in the J1509 tail are in the 
range ~ 10-100 /xG (probably closer to the upper end of this 
range), while the magnetic field in the J1740 tail is a factor of 
a few lower. To obtain more certain estimates, deeper ACIS 
observations and observations outside the X-ray range would 
be most important. 

Using the magnetic field estimates, we can constrain the 
maximum electron Lorentz factor required for producing the 
observed X-ray power-law spectrum: 

7m « 7.6 X 10^(£m/>'m,p)'/^(£/10/uG)-^/2 > 3 x 10*^ (i.e., the 
maximum electron energy Eg^^ > 150 TeV), for the J1509 
tail^. The Larmor radius for such electrons, = 1.7 x 
10'^7(B/10/xG)-' cm, is > 4 x lO"" cm; it is smaller than the 
observed tail width as long as 7^ < 10^". 

The magnetic fields derived above and the measured pro- 
jected tail length allow one to constrain the average projected 
flow speed in the tail: 

Vflow,± ~ l±/Tsyr, ~ 860(/i/l pc)iB/lO fiGf/\E /I keV)^/^ km s"' , 

(8) 

where 

Tsyn = 5.1 X 10S"^B"^ s ~ 1100(5/10 AtG)-^/^(£'/l keV)"^/2 yr 

(9) 

is the synchrotron cooling time [here we used the character- 
istic Lorentz factor, 7 ~ 1.4 x l0^iB/lOiiGT^^\E/lkeVy^^ 
for the electrons that give the main contribution to the syn- 
chrotron radiation at energy E]. Using equation (8), the mag- 
netic fields in Table 4, and the highest energies at which 
the different regions of the J 1509 tail are seen, we obtain 
Vf[ow,_L ^ 7,500 km s*' for region Rl, and Vflow,_L ~ 13,000 
and 12,000 km for regions R1-HR2 and R1+R2H-R3, re- 
spectively, for d = 4 kpc. Given the uncertainty of the mag- 
netic field estimates, the actual average flow speeds can be in 
the range Vflo„._L ^ 5,000-100,000 km s-', with higher val- 
ues more plausible than lower ones. A similar consideration 
for the J1740 tail yields vaow.j. > 1000 km s-', but we should 
note that this is only a conservative lower limit (because the 
actual projected tail length is likely larger than observed), and 
Vfiow = Vflow.i/ sin; can substantially exceed the projected ve- 
locity for small sin;. More certain estimates for vaow could 
be obtained by direct measurements of the proper motion of 
the emission "clumps" (local brightenings) in the J1509 tail 
(see Fig. 5, bottom left), similar to the proper motion mea- 
surements of moving "blobs" in the Vela pulsar jet (Pavlov et 
al. 2003), which would require a series of deep ACIS obser- 
vations. 

Even with account for the large uncertainty, we can con- 
clude that, at least for the J 1509 tail, the flow speed is much 
larger than any reasonable pulsar speed. This means that 

* Sydney University Molonglo Sky Survey 

' We should note that HB07 assumed 7 = 10* and derived the magnetic 
field B ^ 5-7 in the shocked region. Such parameters are inconsistent 
with the X-ray data because they correspond to the characteristic energy of 
synchrotron radiation, E ~ 0.03 eV, in the infrared range (A ~ 30 /im), five 
orders of magnitude below the observed X-ray energies. 
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the equation Tsyn = / / Vpsi , used in many papers on bowshock- 
tail PWNe"', is inapplicable. (This equation impUes that the 
downstream flow instantly acquires the velocity of the ambi- 
ent medium [i.e., vgow = t^'psi- in the pulsar reference frame], 
leaving a synchrotron-emitting "trail" behind the moving pul- 
sar.) On the other hand, the inferred flow speed is smaller than 
I'flow ^ 0.8c-0.9c obtained in numerical simulations by BOS 
for the bulk of the tail volume behind the TS bullet. How- 
ever, those simulation were done only for the beginning of 
the tail (e.g., for distances < l3Ri, from the pulsar in BOS), 
whereas we analyze a much larger extent of the tail (up to 
l± ^ 600/?/,). It seems plausible that the flow is mildly rela- 
tivistic in the immediate vicinity behind the TS bullet, and it 
decelerates at larger distances because of the expansion of the 
overpressurized flow, entrainment of the ISM material (e.g., 
caused by the shear instability at the CD surface separating 
the fast-moving shocked PW from the shocked ambient ISM), 
or internal shocks in the tail. 

The broadening of the J1S09 X-ray tail with increasing 
distance z from the pulsar (Fig. 6) suggests that the cross- 
sectional area of the tail increases up to z « 1 .5', as S(z) oc z'', 
with K ~ 0.7- l.S. At larger distances, 2' <z <^ 5', the tail 
seems to have more or less constant width, although this re- 
quires confirmation with deeper observations. The observed 
radial expansion of the tail at large z significantly exceeds the 
predictions by Romanova et al. (200S; eqs. [32-34]), for plau- 
sible values of vaow and Vpsr (see Fig. 7). Including the inter- 
action between the shocked pulsar wind and shocked ambient 
gas at the CD surface might help to bring the model predic- 
tions in agreement with the observations. Polarimetric obser- 
vations in the radio would allow one to measure the direction 
of the magnetic field in the tail and test whether it is consis- 
tent with the toroidal field assumption adopted in the current 
models We should note that the relationship between B, S, 
and Vflow at given z depends on the magnetic field topology. 
For instance, in ideal MHD models, B cx Vg^^^"'/^ for purely 
toroidal magnetic field (BOS), while B oc if the magnetic 
field is parallel to the tail direction. Deep X-ray observations 
would allow one to measure both S(z) and B{z) and test the 
validity of the ideal MHD approximation and the impact of 
the ambient matter entrainment. 

3.3. Multiwavelength Aspects 

The good alignment between the X-ray tail of J1S09 and the 
elongated radio feature seen in the SUMSS image (Fig. 14) 
suggests their physical association (HB07). As it was noted 
by HB07, the larger length of the J1S09 radio tail could be ex- 
plained by slower synchrotron cooling of the radio-emitting 
electrons. However, the elongated radio feature has highly 
non-uniform surface brightness, with two maxima. The maxi- 
mum located at w 8' from the pulsar is clearly associated with 
an unrelated point source, possibly an AGN. (The source is 
also seen in X-rays; it is marked "R" in Fig. 1.) The origin 
of the other maximum (located about S' from the pulsar) is 
unclear. It appears to be extended and could be a remote SNR 
(as suggested by Whiteoak & Green 1996), a radio galaxy, 
or an inhomogeneity in the pulsar wind flow (e.g., due to an 
internal shock or instability). To discriminate between these 

1° Including HB07. 

" So far, only the Mouse PWN tail has been studied in the radio, and 
the direction of the detected polarization suggests that the magnetic filed is 
predominantly aligned with tail direction (Yusef-Zadeh & Gaensler 2005), 
contrary to the current models. 
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Fig. 14. — Composite X-ray (red) and radio (blue) image of the field of the 
J1509 tail. 

possibilities, deeper radio and X-ray observations are needed. 

Another puzzling aspect is the gradual decrease in the radio 
surface brightness of the tail toward the pulsar (the bright- 
ness drops by at least a factor of S from its maximum value 
at « S' from the pulsar), in contrast with the X-ray surface 
brightness that increases toward the pulsar (see Fig. 14). Un- 
fortunately, the SUMMS image is not deep enough to show 
whether or not the putative radio tail actually connects to the 
pulsar. If the radio emission is indeed associated with the 
J1S09 tail, the decrease in the radio surface brightness to- 
ward the pulsar cannot be explained by the synchrotron cool- 
ing because the corresponding cooling time, r^y*" = 1.7 x 

WiB/lOnGT^IHv/lGnzT^I^ yrs, greatly exceeds the pul- 
sar's age. To explain the lack of radio emission closer to the 
pulsar, one has to assume that the number density of parti- 
cles, n(z), increases with the distance from the pulsar [e.g., 
n{z) oc 1 / Vflow(2) for a constant cross-sectional area of the tail). 
This increase would not result in higher X-ray surface bright- 
ness if the most energetic electrons, responsible for the X-ray 
emission, undergo strong synchrotron cooling, such that their 
emission moves out of the X-ray band by the time when the 
flow reaches the distance where «(z) increases substantially. 
Similar spatial radio-X-ray anti-correlation has been reported 
for the PWNe around the Vela pulsar (Dodson et al. 2003; 
Kargaltsev & Pavlov 2004) and PSR B 1706-44 (Romani et 
al. 200S), which are moving with M.<\. On the other hand, 
in the Mouse PWN (which is the only pulsar tail well studied 
in the radio and X-rays) both the radio and X-ray emission are 
brightest near the pulsar and gradually fade with increasing z- 
Note that the Mouse X-ray tail appears to be shorter (/j^ ^ 1 
pc), or at least its X-ray surface brightness drops faster with 
the distance from the pulsar than in the J1S09 and J1740 tails, 
which could result from faster cooling or lower flow speed. 
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Fig. 15. — Muldwavelength spectra of extended pulsar tails. For J1509 
(top panel), the X-ray and radio fluxes are taken from region R2 (see Fig. 
10). For J1740 {bottom panel; black), the flux is taken from the bright X-ray 
tail region (the ellipse in Fig. 3). For the Mouse PWN (bottom panel; red) the 
X-ray and radio fluxes are for the tail region as defined in the Figs. 3 and 5 of 
Gaensler et al. (2004). An upper limit on the J1740 1.4 GHz flux, estimated 
by V. Kondratiev (priv. comm.) from archival VLA data is also shown in the 
bottom panel. The shaded areas show the uncertainties of the X-ray spectra 
extrapolation. The broken solid lines show examples of possible spectra of 
the J1509 and Mouse tails, which are consistent with both X-ray and radio 
measurements (the detection of the radio tail in J1509 is tentative; see text). 
The spectral breaks correspond to the change of 0.5 in the spectral index due 
to synchrotron cooling. 

We also note that the remarkable 12'-long (Zj^ = 17 pc at c/ = 5 
kpc) radio tail in the Mouse (Yusef-Zadeh & Gaensler 2005) 
is a factor of > 20 longer than its X-ray tail, yet the radio tail 
remains well collimated and straight along its entire visible 
extent. Only deep radio observations with better angular res- 
olution will allow one to understand the nature of the radio 
emission coincident with the J 1509 tail and thus determine 
its true extent. Confirming the radio counterpart of the J 1509 
tail and measuring the magnetic field distribution in the radio 
tail will significantly advance our understanding of collimated 
magnetized outflows from pulsars. 

The multiwavelength (radio through X-rays) spectra of the 
Mouse, J 1509, and J 1740 tails are shown in Figure 15. If 
the tentative interpretation of the extended radio emission 
as a counterpart of the J 1509 X-ray tail is correct, then the 
tail spectrum should exhibit a spectral break between the ra- 
dio and X-ray frequencies. Assuming that the break in the 
spatially integrated spectrum is due to the synchrotron cool- 
ing (Ar ^ 0.5), the break frequency is between ^ 10'^ and 



~ 10'^ Hz, for the allowed range of the X-ray spectral slopes 
(an example of such a break is shown in Fig. 15 for = 2. 16). 
Low-frequency breaks (in a 10-100 GHz range) have been 
previously reported in a number of PWNe, such as G11.2- 
0.3, G16.7H-0.1, and G29.7-0.3 (Bock & Gaensler 2005). The 
spectrum of the Mouse also shows similar behavior, indicat- 
ing a break around 10'^ Hz (Fig. 15). 

No extended radio emission has been reported around 
J 1740. However, this could be due to the fact that the field has 
not yet been observed with sufficiently deep exposures. The 
X-ray spectrum of the J 1740 tail is substantially harder than 
those of the J 1509 and Mouse tails, suggesting a less efficient 
cooling or a different injection spectrum. Also, one should 
keep in mind that, most likely, we have detected only the ~ 2 
pc brighter part of the J 1740 tail, where the cooling makes 
no significant effect, while the entire tail can be substantially 
longer. A deep Chandra ACIS observation is necessary to re- 
veal the true extent of the J 1740 tail and search for spectral 
changes associated with the cooling. 

So far, no TeV emission has been found in the vicinity of 
bowshock-tail X-ray PWNe, with the possible exception of 
the nebula associated with PSR B 1509-58 (Aharonian et al. 
2005), where the origin of the bright feature southeast of the 
pulsar (Gaensler et al. 2002) is not completely understood 
(it could be a jet, a tail, or, more likely, a combination of 
those, i.e., a jet behind the pulsar moving ?A ^ 1). How- 
ever, the X-ray images of J 1509 and J 1740 reveal the pres- 
ence of relativistic electrons with 7 > 10**, which, for a plau- 
sible magnetic field B > 10 /iG, can produce TeV emission 
via the inverse Compton scattering of the background IR and 
CMB photons. Furthermore, given the large angular extent 
of these tails, their TeV emission can be resolved with the 
current TeV observatories such as HESS, MAGIC and VER- 
ITAS. The high energy GeV emission may also be detectable 
with GLAST if there is a substantial density of synchrotron 
radio photons that will be up-scattered into the GeV range. 

3.4. Comparison with other pulsar tails 

Table 5 summarizes the X-ray properties of eight PWNe 
with firmly established long X-ray tails. Two most lumi- 
nous tails belong to the very young, energetic PSR J0537- 
691012 (located in the LMC SNR N157B) and the 26 kyr old 
PSR J 1747-5928 (the Mouse PWN), whose X-ray efficien- 
cies, ?7p„n = Lp„„/E, are 2.2 x 10"^ and 1.6 x 10"^, respec- 
tively. The J 1509 tail ranks third in terms of its efficiency, 
?7pwn ~ 2.6 X lO""*, which exceeds the efficiency of, e.g., the 
Duck PWN (r7pwn = 3.2 x 10"^), created by the younger and 
more powerful PSR B 1757-24, whose spin-down properties 
are very similar to those of PSR J 1747-5928. This once again 
suggests that, although there is a positive correlation of Lpwn 
with E (see Fig. 16), the X-ray PWN luminosity also depends 
on other parameters (see KP08). The remaining four pulsars 
(B1853H-01, J1740H-1000, B0355H-54, and B1929H-10) show 
similar PWN efficiencies, ryp^vn = (2-7) x 10"^, despite the fact 
that their spin-down ages and powers differ by more than 2 
orders of magnitude. Although some of this scatter could be 
attributed to the poorly known distances and different sensi- 
tivities of the individual X-ray observations, it is likely that, 
in addition to the spin-down parameters, there are other pa- 

We should note, however, that, because of high pressure in the young 
SNR, this pulsar may be moving subsonically, in which case the extended 
structure is a "trail" rather than a tail (see §3.2), and hence it cannot be mean- 
ingfully compared with the tails of supersonically moving pulsars. 



X-ray observations of parsec-scale pulsar tails 
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Fig. 16. — X-ray luminosities of PWNe versus pulsar spin-down power 
(see KP08). The data points corresponding to pulsars with firmly established 
long X-ray tails (see Table 5) are shown in color (as well as the names of 
the pulsars). The blue and red points correspond to the total PWN luminosi- 
ties and the luminosities of the compact bright features associated with the 
PWN "heads", respectively. In most cases the luminosity of the extended tail 
component substantially exceeds that of the head. 

rameters governing the radiative efficiency of the winds in 
extended pulsar tails, such as the pulsar speed, the ambient 
pressure, and the angles between the NS magnetic and spin 
axes and the velocity vector 

In addition, for the collimated tails with nearly relativis- 
tic flow speeds, the orientation with respect to the observer 
can be important because the tail brightness is proportional to 
(1+ (3 cos/)'""'"^, where (3 = Vflow/c and / is the angle beween 
the velocity vector and the direction to the observer For in- 
stance, the ratio of the brightnesses at / = 30° (approaching 
flow) and / = 120° (receding flow) can be as large as ~ 10 for 
f3 = 0.3 and T = 2.3. In addition to the brightness changes, the 
outflow components approaching and receding at small an- 
gles will appear shorter simply as a result of projection. This 
might explain why the luminous X-ray tail of the more en- 
ergetic Mouse PWN is so much shorter than that of the less 
energetic J 1509, or why the X-ray tail of the Duck PWN is 
not only short but also very dim compared to the X-ray tail of 
the Mouse. 

Comparing the X-ray luminosities of pulsar tails with those 
of PWNe around slowly moving pulsars, we see that, on av- 
erage, the tails show higher X-ray efficiencies (see Fig. 16). 
This is a natural consequence of the fact that in the ram- 
pressure confined PWNe the energy of the wind flow is chan- 
neled into a narrow cylindrical structure (the tail). This re- 
sults in a higher column density of emitting electrons, hence 
a higher surface brightness, which allows one to detect the 
wind emission at larger distances from the pulsar compared 
to PWNe around slowly moving pulsars where the wind flow 
is more isotropic. 

If the starting flow speed just downstream of the TS back 
surface is more or less the same (mildly relativistic) for dif- 
ferent pulsars, then the observed tail lengths are determined 
by the magnetic field strengths in the tails (because the syn- 
chrotron cooling limits the length) and by the rates at which 
the tail widths increase with increasing distance z from the 



pulsar (because the tail broadening reduces its surface bright- 
ness). The magnetic field in the tail should be proportional 
to (aEy/-^ (a is magnetization parameter in the pre-shock 
wind). The lack of correlation with E in the current, very 
limited sample (Table 5) might be attributed to different wind 
magnetization in different pulsars and/or a large spread of the 
unknown inclination angles of the tails. The tail width is ex- 
pected to show a negative correlation with the ambient pres- 
sure (e.g., /?taii oc [E/ivl^^p^^b)]^/^z^/^ in the Romanova et 
al. 2005 model); however, the dependence is too weak to be 
tested with the existing data. 

4. SUMMARY 

We have presented a detailed analysis of the properties of 
two PWNe associated with the supersonically moving pulsars 
J1509-5850 and J1740H-1000. In both cases, the X-ray im- 
ages reveal parsec-scale tails connected to the pulsars. The 
observed extraordinary lengths of the tails suggest that the 
tail flow starts as mildly relativistic, and its speed remains 
much higher than the pulsar speed up to a distance of a few 
parsecs. The large extent of the J1509 tail suggests an aver- 
age flow speed in excess of ~ 5000t/4 km s"' (a conservative 
lower limit). We estimate the average equipartition field in 
the J 1509 tail to be a few xlO~^ G, with a hint of weaken- 
ing with the distance from the pulsar. The equipartition field 
in the J 1740 tail is likely to be a factor of a few lower A 
noticeable difference between the two tails is that the spec- 
trum of the J 1740 tail is substantially harder, implying that 
either the cooling is not as efficient in the J 1740 tail (con- 
sistent with the inferred lower magnetic field) or the electron 
injection spectrum of J1740 is intrinsically harder compared 
to that of J 1509. The projected orientation of the J 1740 tail on 
the sky suggests that the pulsar is moving toward the Galac- 
tic plane. This means that the pulsar was born well above the 
Galactic disk rather than ejected from the disk with a high 
speed. 

The existing shallow radio survey data suggest a possible 
radio counterpart to the J1509 tail; however, a confusion with 
the background sources remains a possibility until better qual- 
ity radio data are obtained. If the radio feature is indeed as- 
sociated with the J 1509 tail, then the tail spectrum should ex- 
hibit a break between lO'^ and 10'^ Hz. 

The comparison between the J 1509 and J 1740 tails and the 
X-ray tails from other fast-moving pulsars shows that the X- 
ray luminosity does not correlate well with the pulsar spin- 
down luminosity and age; other physical parameters such as 
the pulsar velocity, the ambient pressure and sound speed, and 
the angle between the pulsar's spin and magnetic axis may be 
equally important. In general, we find that the X-ray efficien- 
cies of the ram-pressure confined PWNe are systematically 
higher than those of PWNe around slowly moving pulsars 
with similar spindown parameters. 
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TABLE 5 
Properties of long pulsar tails. 



PSR 


SNR 




vx" 


logr 


logE 




logLpwH** 


log»?pwii 


Rad.^ 


Ref.f 






kpc 


km s~' 


[yrs] 


[ergs s"'] 


pc 


[ergs s"' ] 








J0537-6910 


N157B 


50 


600 


3.70 


38.68 


3.7 


36.43 ± 0.02 


-2.24 


Y 


1,2 


B 1757-24 




5 


130 


4.19 


36.41 


0.5 


33.20 ±0.14 


-3.21 


Y 


3 


J1747-2958 




5 


300 


4.41 


36.40 


1.1 


34.75±0.15 


-1.65 


Y 


4 


J1509-5850 




4 


400 


5.19 


35.71 


6.5 


33.04 ±0.05 


-2.67 


P 


tw 


B1853+01 


W44 


3 


100 


4.31 


35.63 


1.3 


32.39±0.10 


-3.24 


Y 


5 


J1740+1000 




1.4 


400 


5.06 


35.36 


2 


31.54 ±0.20 


-3.82 


N 


tw 


B0355+54 




1.04P 


61^12 


5.75 


34.66 


1.5 


31.49±0.10 


-3.16 


N 


6 


B1929+10 




0.36P 


177!^ 


6.49 


33.59 


1.5 


29.87 ±0.25 


-3.72 


P 


7,8 



Note. — The table includes pulsars with established parsec-scale X-ray tails. 

" Adopted distance to the pulsar (superscript P marks the distances obtained from parallax measurements). 

''Pulsar's transverse velocity, either obtained from proper motion and parallax measurements (for B0355+54 and B 1929+10; Chatterjee et al. 2004) or estimated 
from eq. (3), with n = 1 cm""* and R/, measured from the X-ray data. The exception is J0537-6910, for which we quote the velocity estimated by Wang et al. 
(2001) from the pulsar travel time (this pulsar may be moving subsonically because of the high sound speed in the young host SNR). Also note that since PSR 
B 1853+01 is embedded in a plasma that shows strong thermal X-ray emission, the assumption n = 1 cm"' likely overestimates vx in this case. On the contrary, 
since PSR J1740+1000 is likely moving in a rarefied medium above the Galactic plane, its speed is probably underestimated (see text for discussion). 

"^Projected length of the tail. 

''X-ray luminosity of the tail (including the brighter compact component in the pulsar vicinity) in the 0.5-8 keV band as we measured it from Chandra and 
XMM-Newton data. The values may differ from those published in the original works (see the references in the last column). 

'^Is the PWN detected in radio? P = possibly. 

'References to the original papers where the corresponding Chandra and XMM-Newton data have been analyzed. - (1) Wang et al. 2001; (2) Chen et al. 2006; 
(3) Kaspi et al. 2001; (4) Gaensler et al. 2004; (5) Petre et al. 2002; (6) McGowan et al. 2006; (7) Becker et al. 2006; (8) Misanovic et al. 2007; tw = this work. 
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